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ABSTRACT 

We propose a new electromagnetic-emission mechanism in magnetized, force-free plasma, which is 
driven by the evolution of the underlying dynamic spacetime. In particular, the emission power and 
angular distribution of the emitted fast-magnetosonic and Alfven waves are separately determined. 
Previous numerical simulations of binary black hole mergers occurring within magnetized plasma 
have recorded copious amounts of electromagnetic radiation that, in addition to collimated jets, in¬ 
clude an unexplained, isotropic component which becomes dominant close to merger. This raises 
the possibility of multimessenger gravitational-wave and electromagnetic observations on binary black 
hole systems. The mechanism proposed here provides a candidate analytical characterization of the 
numerical results, and when combined with previously understood mechanisms such as the Blandford- 
Znajek process and kinetic-motion-driven radiation, allows us to construct a classification of different 
electromagnetic radiation components seen in the inspiral stage of compact-binary coalescences. 
Subject headings: gravitation — gravitational waves — plasmas — radiation mechanisms: general 


1. INTRODUCTION 

With the imminent direct detection of gravi tational 
waves (GW s) by second generation detectors (Dooley 
et al.|(|2015|)), the pursuit of an understanding of the elec- 
tromagnetic (EM) counterparts to GWs becomes urgent, 
as a joint observation in both channels will provide irre¬ 
placeable m eans to diagnose propertie s of the astrophys¬ 


ical sources (Christensen et al. 
important types oi sources tha 


( 2011 )). 

could r 


One of the most 
adiate both grav¬ 


itationally and electromagnetically is a coalescing com¬ 
pact binary, involving black holes and/or neutron stars 
surrounded by magnetized plasma (forming the so-called 
“magnetospheres”). The magnetic field could originate 
from the accretion disk of the binary or neutron stars 
themselves, and the plasma could be generated from vac¬ 
uum polarization, and/or charged particles coming off of 
the star surfaces and the accretion disk. Recent numeri- 


cal simulations ( 

3 alenzuela et al. 

(2010b); 

Neilsen et al. 

(2011) 

Alic et al 

(2012)) have shown that EM radiation 

is indeed given c 

ft by such systems in abundance even 


before merger and for binary black hole systems (while 
current joint-observation efforts concentrate on the post- 
rnerger stage of system s with at least one neutron star 
(Nissanke et al. (20131)), providing further optimism for 
the success of multi-messenger astronomy. The next step 
is then to clarify the various physical processes at work 
that, together, produces the EM signals seen numerically 
(in particular, an isotropic radiation that dominates near 
merger time has not been previously understood analyt¬ 
ically). A complete classification and characterization of 
these processes is a prerequisite for extracting useful in¬ 
formation about the binary systems from the observed 
EM signals. We provide such an analytical clraracteriza- 
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tion in this work and compare it with previous numerical 
results (see Fig. [I] below). 

Within magnetospheres, the energy density of the mag¬ 
netic field often dominates over that of the plasma parti¬ 
cles, creating what’s referred to as a force-free plasma. 
Thank s to the seminal works b y Goldreich fe Julian| 
(1969) and Blandford & Znajek (1977p, it is widely ac- 
cepted that iorce-free plasma can act as a medium for 
powering outgoing EM radiation (or jets) at a cost of re- 


holes 

Thorne (1994) 

Meier 

(2012); Palenzueia et al. 

(2011) 

; Spruit et afo( 

997); Hansen & Lyutikov 

2001)). 

iVlore recent studies (Hansen 

& Lyutikov (2001 

Lyu- 

tikov 

2011): Brennan & Gratia (2013): Henna 

(2015); 

Palenzueia et al. [2009, 2010c b a); JNeiisen et al. 

(2011); 

iVloesta et al. (2012); Alic et al. (2012); Faschaiid 

is et al. 

(2013) 

(2013) 

; Me Wiliams & Levin 

(2011); DUrazio & Levin 

; [Morozova et al. (2014 

)) suggest that a force-tree 


ergy of moving objects to power EM radiations in the 
form of jets launching from star surfaces (or the black 
hole horizon), accompanied by some isotropic flux. We 
refer to this as the kinetic-motion-driven radiation [^] 
which i s also seen from satell ites moving in earth’s iono¬ 
sphere (fDrell et al. (1965a|b|). 

There is however/a third mechanism, which we shall 
call the gravitation-driven radiation , which will be the 
focus of this paper. When the background spacetime 
becomes dynamic, the local EM energy density of mag¬ 
netized plasma deviates from its equilibrium values and 
these inhomogeneities tend to propagate out via plasma 
waves. A similar phenomenon is known to exist in space- 
times without matter (the Gertsenshtei n-Zeldovich effect 
(Gertsenshtein (1962); Zeldovich (19731)), where the out- 


5 We caution that although terms like this have been used, here 
and in other literature, to label different EM emission mechanisms, 
the nonlinearity of force-free dynamics and gravity makes a math¬ 
ematically rigorous classification difficult. This is particularly true 
with fully nonlinear numerical simulations, in which all of the emis¬ 
sion mechanisms discussed here and below are likely present, even 
when the simulation aims to study a particular one. 
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going radiation consists purely of vacuum EM waves. 
In addition, the generation of magnetohydrodynamic 
(MHD) waves by the influence of gr avitational waves has 
been examined in Duez et al. (2005). Although this effect 
has not been explicitly discussed in the context of force- 
free magnetospheres, we note that force-free electrody¬ 
namics (FFE) can be viewed as the l ow-inertia limit of 
relativistic magnetohydrody namics (McKinney (2006); 


Paschalidis & Shapiro (2013)). In this paper, we will ex- 
amine essentially the same physical process, but where 
the driving gravitational dynamics is not an (idealized 
wave-zone) gravitational wave. Within force-free plasma, 
energy can be carried away by two different classes of 
waves. One class is called the fast-magnetosonic waves in 
the local short-wavelength limit (the wavelength is much 
smaller than the radius of spacetime curvature), whose 
global and longer-wa velength counterparts a r e named the 


“trapp ed modes” in Yang & Zhang (2014); Yang et al. 
(2015). These tend to behave similarly to vacuum EM 
waves and propagate in a more egalitarian fashion in 
terms of sky directions. The other class of waves a re th e 


Alfven waves , generalizing to “tra veling waves” (Yang & 
Zhang (2014|); Yang et al. d2015|S or principal nu ll so¬ 


lutions (|h5renna^^fari2013|); Zhang et al. (2015)). A 
salient feature of the Alfven waves and their generaliza¬ 
tions (for brevity, we will not distinguish between them 
below, similarly for the other class) is that they propa¬ 
gate along the magnetic field lines, and are as such auto¬ 
matically collimated if the magnetic field threads through 
the orbital plane of the binary nearly orthogonally (a 
natural configuration for accretion-disk-supported field). 
Below, we show how to compute their fluxes as generated 
by the gravitationally-driven process. 

In order to perform the a nalysis, we apply the geo ¬ 


metric approach promoted by Gralla &; Jacobson 
Uchidal (fl997a] |l998[ |1997b| )r lHarter| dl979| );' 


(20141; 


E chida 


( 1997c|d I, whose introduction has t riggered many new 


develo p ments ([Zh ang et ah ([ 2014 ); Gralla fe Jacobson 


; Lupsasca ct al. (2014); Lupsasca &; Rodriguez 


; Gralla et al.| (]20i5a|b|)) m obtaining exact solu¬ 
tions to force-tree electrodynamics, in addition to new 


interpretat i ons of previous results (Penna (2014); Menon 


& Dermer (2007)). Unless otherwise specified, the for- 


mufae below are in natural units, so that c = G = 1. 


2. SET-UP OF THE PROBLEM 


Let us assume that there is a stationary FFE config¬ 
uration in a stationary background s paceti me with met ¬ 
ric (?b. According to dis cussions in |Uchida| (|1997c|d|); 
Gralla & Jacobson (2014), it is possible to find at least 
one pair oK“Euler potentials” <(> 15 , 25 , such that F-q = 
dip ib A d(f> 2 B, where Eg is the background electromag¬ 
netic field tensor. Now suppose that the spacetime be¬ 
comes dynamic and its metric is g = <?b + £0, where e 
parametrizes the magnitude of the spacetime deforma¬ 
tion from its stationary state. Correspondingly the Eu¬ 
ler potentials will also deviate from their original values: 
01.2 = 0i,2B + £001,2, wh ereby the non-linear FFE w ave 
equations they satisfy are Gralla & Jacobson (2014) 

( 1 ) 


d0i,2 A d* F = 0 , 


with F = 00i A 00 2 . Note that the Hodge star * is 
now with respect to the total metric g , so that it de¬ 
pends on metric perturbations. In order to study the 


gravitationally-induced plasma waves, we shall linearise 
the above equation to the leading order in e, and obtain 


0001,2 A0*b Eb+00i,2B A 0 *b SF 

d * Fa 

= — 001,2B A 0——- . (2) 

at 

This equation describes the excitation of plasma fields 
00 i ,2 by the source on the right hand side, which is linear 
in h. It implies that GWs interacting with magnetized 
plasma can generate plasma waves. Moreover, it predicts 
that a time-dependent Newtonian source within magne¬ 
tized plasma also induces plasma radiation, an effect that 
has been overlooked before and could have observational 
consequences. 

3. RADIATION IN NEARLY FLAT SPACETIMES 

Now we specialize to a simple yet important exam¬ 
ple where the background metric is flat, i.e., g^ = 
Vfj.u + This is a good approximation when the 

gravitational field generated by matter sources or GWs 
is weak. In addition, let us assume that the plasma is 
magnetized along the z direction, with field strength B 
so that Fb = Bdx A dy. When the spacetime becomes 
dynamic, the EM field 2-form can be written as (note 
we consider only those FFE perturbations driven by the 
spacetime variations, and so use the same flag e) 

F = B(dx + e00i) A (dy + eS(f> 2 ). (3) 

With this set-up, one can straightforwardly work out the 
Hodge star rules, plug them into Eq. ([2 ), and obtain a 
coupled set of wave equations for 00i,2- These equations 
can further be diagonalized through the definition of a 
new set of variables: 


01 = <9*002 - <9y001, 02 = <9y002 + <9*001 , (4) 

in which case the wave equations decouple into 

(-B 2 + - 32 htx _ 32ht y 4 - 32h y* - 32flxz 

[ 4 dtdy dtdx dzdx dzdy ’ 

( —+ V")02 = -(<9* + + hyy — h zz ) 

d 2 h yz d 2 h xz d 2 h yt d 2 h xt 
dydz dxdz dydt dxdt 

The first equation describes a wave propagating along the 
magnetic field lines, or in other words the Alfven wave. 
The second equation describes the fast-magnetosonic 
wave, which propagates in all directions. These equations 
are gauge-invariant, as can be checked by substituting in 
the infinitesimal gauge transformation Xi —> Xi + that 
leads to 

hyii y hyv -)- (0 ji/ -(- ^u\/l ~ “t“ (0,/X , (6) 

and 

001^ 001 +£*, 002 ^ 002 +^- (7) 

Denote the source terms in Eq. [5] as Si and 5 2 respec¬ 
tively, the solutions to these wave equations can be ob- 
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tained through the use of Green’s functions, 


4>i = 


1 

2 


dzdt Q(t — t! — \z — z'\)Si(t', z ’), 




j3„ S 2(t~ |X-X'|,X') 

47t|x — x'| 


( 8 ) 


where 0 denotes the Heaviside step function. After eval¬ 
uating ipi t 2 , we can reconstruct 5(f> and subsequently 
F , by noting that 


(dl + d 2 )S(p 2 = d x ip! + dyip 2 , 

{dl + = -dyip! -(- d x ip 2 , (9) 

whose solutions are (applying the Green’s function for 
2-D elliptic equations) 


x , [ a >a ( x ',y',z) + A y ip 2 (x',y\z) 

5h = J dxdy -MiJTip- 

/A..' A ^‘ 2 ( x ',y',z) - A y ipi(x',y',z) 


S(j )i = / dx dy‘ 


2ti-(A2 + AS) 


, ( 10 ) 


4. FLUX EXTRACTION 

According to Eq. (p|, the fast-magnetosonic waves are 
quite similar to the vacuum EM waves, where the source 
term S 2 can also be decomposed into multipole contribu¬ 
tions. Let us assume that the binary (with total mass M) 
is practicing near-circular motion with a period of 2i r/f2, 
in which case ip 2 in the radiative zone can be written as 


V>2 


£/*,(«>) 


gim[<f>— Q(i—r)] 

r 


( 12 ) 


where the m = 0 piece corresponds to the DC monopole 
field, which does not radiate. The coefficients f m may 
be further decomposed into a summation of associated 
Lengendre polynomials, starting from l > \m\. In order 
to compute the energy flux, we need to reconstruct 8<p 1,2 
with Eq. ( fl0| (in the absence of ip\), or more efficiently, 
by noticing t hat 8(f) 1,2 must possess similar asymptotic 
forms as Eq. (12): 


5(f>ia ~ X^m 2 (0) 


0 im[(p—Q(t—r)] 


(13) 


where = x — x' and A y = y — y'. 

Analogous to the Gertsenshtein-Zeldovich effect, 
Eq. ([5]) together with Eq. ([8]) explicitly show that GWs 
injected into magnetized plasma would generate both 
Alfven and fast-magnetosonic waves. Suppose that the 
gravitational wave packet has a characteristic amplitude 
h and a length-scale of A, it is then straightforward 
to see that the plasma-wave luminosity £qw satisfies 
£gw cx B 2 \ 2 h 2 \ a relationship that can be compared 
with future numerical experiments. Here we focus in¬ 
stead on the case where the source is generated by two 
orbiting compact masses, in order to study the radiation 
of a binary system in the inspiral stage. With a Newto¬ 
nian matter source (as the leading order post-Newtonian 
term of general relativistic expre ssions, which i s suffi cient 
for our purpose), h is given by (Misner et al. (1973)) 


hnn — 2 d’ 


3 / P( x ') 


hjk — 28 


jk 


d x 


X - x'l 

/ p(x') 


h 0j = -4 Id' 


3 , p( x 'K 


( 11 ) 


When the source consists of a pair of orbiting black 
holes, the formulae above are valid at places away from 
the black holes, which are themselves replaced by point 
masses. However, the Newtonian approximation be¬ 
comes inaccurate near the black holes. In addition, in 
order to compute the plasma waves at far away and ex¬ 
tract the energy flux, we must exclude the points en¬ 
closed by the black hole horizons. Therefore, in practise 
(see Sec. [5]), we remove two excision spheres when com¬ 
puting the integrals in Eq. [8j To test the sensitivity of 
the gravitation-driven luminosity values on the excision 
radii choice, we vary their values from 2 M to 3 M (M be¬ 
ing the black-hole mass), and observed that the resulting 
flux changes less than 10 precent. For the presentation 
of data in Sec. [5] then, we adopt the cut-off radius choice 
of 3 M. We caution that this insensitivity to excision 
radii could change significantly if we take into account 
relativistic (Post-Newtonian) corrections to the metric. 


and the relationship between gf^ 2 and f m can be obtained 
using Eq. Q with if) 1 = 0: 

9m\&) = ~ n — 0 • fl[/m+lW + fm-l(d)] , 

2mi l sin 0 


We can then substitute these expressions into 8<f> 1 and 
8(f> 2 , and subsequently Eq. pi to obtain the field 2-form. 
It is then straightforward, although tedious, to extract 
from it the electric and magnetic field vectors, and com¬ 
pute the Poynting vector. In the end, we arrive at the 
flux formula for fast-magnetosonic waves: 


(9) = - 


1 


2mil sin 9 


[f m+ i(e)-f m -i(e)\. (14) 


*$fast — 


B 2 \fm(0)\ 2 CSC 2 (6) 


mzfi 0 


(15) 


The Alfven waves, on the other hand, propagate along 
the magnetic field lines. Based on Eq. pi), we write if) 1 
in the radiative zone \z\ M as 


V’l = y dk x dkyA ± (, k x , ky, u, i))e ikxX+ikyV . (16) 


where ± stands for the top/down extraction surfaces and 
u= t — z, v = t + z. The effective radiative part of ip 1 
is only a function of u for z M, and a function of v 
for M. One can write the associated Scpi^ in a 

similar format, which satisfies Eq. with ^2 = 0: 

to- J dk * dk >j^Tki A±e “" x * lt "‘• 

A ± = Af(k x ,k y ,v) + A±(k x ,k v ,u ), (17) 


from which we obtain the luminosity function 

±\d u A±\ 2 T \d v A±\ 2 


Ca if = 2 B 2 Y, 


dk x dk y - 


-X- z-2 
^x i ^y 


(18) 
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For systems with mirror symmetry about the orbital 
plane, it suffices to only compute the luminosity on one 
side and double the result. 

5. BINARY BLACK HOLE COALESCENCE 

We can now compare our analytical predictions with 
numerical simulations of equal-mass binary black hole 
coalescences, and try to identify the physical mecha¬ 
nisms behind the “isotropic” and “collimated” EM radi- 


ations seen there (iPalenzuela et al. (2010b I; 

Nfeilsen et al. 

(2011) 

Alic et al. 

(2012); Brennan & Grail 

1 (2013)), as 


itate comparison, we adopt the same contextual parame¬ 
ters as in the numerical experiments above, i.e., a binary 
black hole system with 10 8 solar masses for each hole and 
a background magnetic field at 10 4 Gauss. We also note 
that the strength of the EM emissions is much weaker 
than that of the gravitational-wave emission, where the 
gravitational radiation-reaction leads to the shrinking of 
the orbital radius. As a result, it is a valid and com¬ 
mon approximation to ignore any back-reaction of the 
EM radiations on the evolution of the spacetime. 

Both fast-magnetosonic and Alfven waves are pro¬ 
duced during the sequence (inspiral, merger, and then 
ringdown) of binary merger stages, and they radiate 
mostly in the forms of “isotropic” and “collimated” 
fluxes, respectively. Below, we will concentrate on the 
inspiral stage (leading into the merger itself) that’s the 
most interesting for multi-messenger astronomy. Dur¬ 
ing this stage, the EM emissions can be classified into 
rotation-driven, kinetic-motion-driven and gravitation- 
driven types. The rotation-driven radiation is generated 
by the Blandford-Znajek mechanism, which sup ports a 
jet- like radiation with luminosity of the order (Neilsen 

eFahlpOTII)) 


C r 


2.4 X 10‘Ws (j|c) (loSj) < 19 > 

in cgs units and when spin is aligned with the magnetic 
field, or abbreviated as 2.4L 43 Bf M)|d 2 . Here Mi is the 
ith black hole mass and eh is the dimensionless spin pa¬ 
rameter of the black hole ranging from 0 to 1. 

As a black hole moves through magnetized force- 
free plasma, it l aunches collimated jets a l ong the mag¬ 
netic f ield lines (Palenzuela et al. (2010b); Neilsen et al. 

(20111). The power of this radiation is proportional to v 2 
and thus H 2 / 3 . In addition, if the black hole also follows 
accelerated motion, it generates an additional Poynting 
flux similar to accelerated charges in vacuum, which can 
be attributed to fast-magnetosonic wave emission. Its 
power is on the order of 2/3 q 2 a 2 (“Larmor formula” of 
Brennan & Gralla (2013)), where the effective monopolc 
charge q should have value ~ 2 BM 2 and the acceleration 
obeys a oc v 2 /d oc fl 4 / 3 . Summing up the two contribu¬ 
tions, we have for kinetic-motion-driven radiation that 

£ m - 1 .6L 43 B 2 Ml /3 n 2 J* + 0.5L 4 o-BfMg 4/3 S2^4 . (20) 

The merger happens at around S2 ~ 2 x 10 _4 s _1 , and 
so the acceleration-induced radiation is sub-dominant 
through the entire inspiral stage. 

We now turn to t he g ravitation-driven radiation. With 
Eqs. (§, (§. and ( |ll[ ), we can estimate the orbital fre- 



Log 10 (n) 

Fig. 1.— The total fast-magnetosonic £f ast (blue, in units of 
ergs/s), Alfven wave Cau (purple), an d motion-driven (orange) 
luminosities as functions of the orbital frequency £1 (in units of 
1/s). For £fast and ^Aif 5 the dots represent numerical integration 
results; the dashed lines are linear fits with flexible slopes; the solid 
lines are fits with a fixed slope of 4/3. For £ m , the curve is from 
Eq. ( |20| ) (although both terms are included, only the first term 
corresponding to collimated radiation is significant). For numerical 
data, the c rosses are the meas ured isotropic and collimated fluxes 
taken from Neilse n et al.| ( |201 11) for the non-spinning binary system. 

quency dependence of this class of EM emissions for 
a binary black hole system. The source term of fast- 
magnetosonic waves scales as M/d 3 , where d is the 
orbital separation. Such a source term ge nera tes -02 
in the multipolar-expansion manner of Eq. (12), with 


the luminosity for each multipole moment scaling as 
B 2 M 2 v 21 oc fl 2Z / 3 . For unequal mass binaries, the ra¬ 
diation contains a dipole piece with 1 = 1, whereas emis¬ 
sion from an equal-mass binary starts at the quadrupolar 
order (l = 2). On the other hand, the source term for 
Alfven waves scales as Mvil/d 2 and the corresponding 
flux scales as B 2 M 2 v A oc H 4 / 3 . 

In Fig. [I] we plot the fl-dependent luminosities for 
both fast-magnetosonic and Alfven waves, foiyan equal- 
mass b inary system (as is sim ulate d in |Neilsen et a l. 
(20111, Palenzuela et al. (2010bI, and Alic et al. (12012 (J) 
with the cut-otl radius chosen at 1.5 times the horizon 
radius (it turns out that the results are insensitive to the 
cut-off radius). More specifically, we substitute the den¬ 
sity profiles appropriate for poin t m asses following New¬ 
tonian Keplerian orbits into Eq. (11), and feed the result¬ 
ing metric perturbation into the light hand side of Eq. ([5j 
to obtain the expressions for Si and S/• These then allow 
us to numerically integrate out Eq. ([8]) and acquire ip i 
and ip 2 , representing the Alfven and fast-magnetosonic 
waves respectively. To compute the Alfven flux Haie we 
apply d u and d v to ipi and take the results through a 
numerically Fourier transformatio n pr ocedure to obtain 
d u A ± and d v A ± according to Eq. (16). Finally, another 
numerical integration according toTlq. (18) provides us 
with £aie We do this for several black hole separations, 
as signified by their different Keplerian orbital frequen¬ 
cies, and plot the results as the purple dots in Fig. [T| We 
also compute the fast-magnetosonic fluxes £f as t at these 
separations. In this case, we simply need to project rip 2 
onto exp (imp) basis (taking m up to 30) and substitute 
the resulting f m values into Eq.|l5]to compute £f as t- The 
results are plotted as the blue dots in Fig. [l] 

From the figure, we can see that the luminosity values 
are consistent with the quadrupolar contribution’s dom- 
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inance over higher multipoles, with a fl 4 / 3 scaling. We 
can also read off the dependence of Cam and £f as t on 
the m agne tic he ld s trength from their respective formula 
(Eqs. (181 and ©), which is .Bf. Simple dimensional 
consideration fixes the dependence on Mg for us, which 

is Mo 10/3 . What remains to obtain a formula similar to 
Eq. (201 for the gravitation-driven case is the determina¬ 
tion of the coefficients of proportionality, which set the 
overall amplitudes for the fluxes. These are simply the 
intercepts on the vertical axis of the solid purple and 
blue fitting lines in Fig. |T] (in other words, they come 
from actually solving the equations and are not new in¬ 
dependent rough estimates). In the end, we obtain that 
the gravitation-driven radiation should scale as 


Cq —-Tfast T Cam 

wl.7. L 43 BlM^ 0/3 fl^ + 0.58L 43 B|M 8 10/3 fll / 4 3 ■ 

( 21 ) 

Close to merger, the gravitation-driven, fast magne- 
tosonic radiation dominates over flux contributions from 



account nonlinearities, so the analytical fit to numerical 
data should be interpreted with a pinch of salt. The aim 
of the present paper is only to demonstrate the existence 
of the gravitation-driven radiation, and the fact it can po¬ 
tentially produce large fluxes, especially an isotropic one 
during merger, rather than trying to make a fit to the nu¬ 
merical data with our zeroth-order calculation. In partic¬ 
ular, our results should in no way be interpreted as fully 
“explaining” the numerical results. In particular, we note 
that the matching for the fast-magnetosonic/collimated 
flux (blue crosses versus blue lines) at low frequencies is 
less accurate. Without a detailed examination involving 
targeted numerical experiments and higher order ana¬ 
lytical computations, we can not state with certainty the 
exact reason for this, so future studies are required. Here, 
we can but point out some more obvious subtleties in the 
matching procedure. 

Most importantly, as mentioned above, the Newtonian 
approximation breaks down in the vicinity of the black- 
holes in our zeroth-order calculation, and this happens 
regardless of the orbital separation. Although the fluxes 
change by only a few percentage points when we move 
the inner cut-off radius from 3M to 2M, the dominant 
contribution to our numerical flux integrations neverthe¬ 
less originate from the neighbourhoods of the black holes, 
instead of the wavezone. Therefore the omission of non¬ 
linear relativistic effects might be the main approxima¬ 
tion here, and taking into account the post-Newtonian or 
relativistic corrections may further change the luminos¬ 
ity estimates above. Other effects, such as the absorption 
by black holes, should also be treated properly. 

Secondly, the numerical fluxes are divided according 
to their directions of propagation, catering more for the 
observational consequences than for matching with an¬ 
alytical classifications. Such imperfect correspondences 
between concepts employed by numerical and analytical 
studies lead to systematic matching errors. For exam¬ 


ple, the collimation in the numerical study is defined to 
be flux propagating inside a cone of a certain opening 
angle, in analogy with the usual jet language, while for 
Alfven waves climbing the vertical magnetic field lines, 
a cylinder enclosing the binary (or two cylinders around 
individual black holes when they are far apart) would 
be more appropriate. Therefore, with a large extrac¬ 
tion radius and when the black hole separation is large, 
the numerical cone would likely enclose a fair amount 
of fast-magnetosonic waves, contributing to the relative 
weakness of numerically measured isotropic flux. Many 
other numerical difficulties associated with subtracting 
off a background radiation in order to construct a divi¬ 
sion of the overall flux into the collimated and isotropic 
types, especially when the overall flux is weak , have been 
discus sed in the numerical paper s such as |Neilsen et al.| 
(20111 and Moesta et al. (2012|). We refer interested 
readers to these important literature. 

In the future, more specifically designed numerical ex¬ 
periments are necessary to test this gravitation-driven 
emission mechanism, including possibly binary star, in¬ 
stead of binary black hole, simulations. Improved so¬ 
phistication in analytical computations is also necessary, 
before the effects of the various simplifying assumptions 
we made in the present work can be disentangled. 


6. DISCUSSION 

We briefly comment on plasma wave generation during 
the other stages of binary black hole coalescences. Dur¬ 
ing the merger phase, both the spacetime and the magne¬ 
tosphere are highly dynamic, and the best tool to under¬ 
stand their evolution is through numerical simulations. 
However, in the ringdown stage, the time-dependent part 
of the emission arises from: (i) the ringdown of the 
magnetos p here, as described by its eigenwaves (Yang 
& Zhang| (|2014|); [Yang et al.l (|2015|)); (ii) the gravita- 
tional quasmormal modes will drive additional emission 
by coupling to the stationary part of the black-hole jets, 
an effect quantifiable using black-hole perturbation the¬ 
ory. Note that by the “ringdown” stage, we mean the 
period before the post-merger black hole settles down 
to Kerr. The settling time can be estimated as l/w| 2 , 
where w| 2 is the imaginary part of the frequency for the 
l = 2, to = 2 quasinormal mode (the dominant mode). 
The value of w 22 is about 0.1/M for Scliwarzschild black 
holes and w| 2 ~ \/l — a/M for rapidly spinning black 
holes, which asymptotes to zero in the extremal spin 
limit (i.e. the modes are long lived and the settling is 
protracted) r] For a post-merger black hole of 10 s so¬ 
lar masses, the Scliwarzschild formula translates into a 
settling time of about eight and a half hours. So al¬ 
though extremely transient in nature, this period may 
be observationally detectable. On the other hand, the 
real part is ~ 1/M for rapidly spinning black holes 
and ~ 0.5/M for Scliwarzschild black holes. During the 
ringdown stage, the gravitation-driven luminosity can be 
estimated as 


C G ~ 5 2 M 2 (c^ 2 M) 4/3 e~ 2 ^2* ? (22) 


6 Fo r generic Kerr black holes, please see Fig. 5 in |Yang et al.| 
( |2012[ > for the mode decay rates. 
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while the Blandford-Znajek flux is approximately 

£bz ~ B 2 M 2 (af/M) 2 , (23) 

where a,f is the spin parameter for the final black hole. 
As the final black hole in generic binary mergers is rotat¬ 
ing, we expect the Blandford-Znajek contribution to be 
important, and the gravitation-driven emission to also 
be an important part of the total flux at least within a 
timescale of l/w| 2 - 

During the ringdown stage, both the spacetime met¬ 
ric and the magnetosphere would be time-dependent, 
with simi lar but not exactly th e same characteristic fre¬ 
quencies (Yang & Zhang (2014)). The gravitation-driven 
mechanism would account tor the metric variation’s mod¬ 
ifying effect to e.g., the Blandford-Znajek process, but 
not that from the magnetosphere ringing. In other words, 
multiple transient effects are present and it would be dif¬ 
ficult to disentangle the signals they generate. Never¬ 
theless, if quasi-periodic flux variations from the post¬ 
merger black hole can be detected, then one could in 
principal do interesting measurements such as that on 
the black hole spin. 

For completeness, we can also estimate the flux mod¬ 
ification due to the presence of current-sheets near the 
black holes, which is approximately the geometric mean 
of collim ated and acceleration-in duced radiations (see 
Eq. 42 in Brennan & Gralla (2013)). With units restored 


and according to Eq. (20), the corresponding luminosity 
is sub-dominant near merger. In addition, although we 
have examined the gravitation-driven plasma wave gen¬ 
eration here in the context of force-free plasma, we expect 
similar signatures to persist in materials following more 
generic MHD equations. 

Finally, we note that in the binary black hole exam¬ 
ple, energy is emitted at very low frequencies (below the 
plasma frequency). In fact, during the Blandford-Znajek 
process, the outgoing energy flux is carried out at the 
DC frequency. This is allowed for MHD waves (includ¬ 
ing waves in f orce-free plas ma), but not for unmagnetized 
plasma (Thorne & Blandford (2016)). 
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